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Abstract 

We present preliminary results from a search for inclusive charmless B — > KX decays. These 
decays occur dominantly via one-loop b — > s penguin transitions, and can provide useful information 
about these processes. Using a sample of 288.5 fb _1 collected with the BABAR detector at the 
PEP-II asymmetric-energy e + e~ B Factory at SLAC, we search for high-energy kaons recoiling 
against fully reconstructed B decays. We measure the partial branching fractions for kaons with 
momentum p*(K) > 2.34 GeV in the B rest frame, and obtain (in units of 10 -6 ): B(B — ► K + X, p* > 
2.34GeV) = 196i^(stat.)±^(syst.) and B(B -> K°X,p* > 2.34GeV) = 154lg(stat.)l^(syst.) 
(< 266 at 90% C.L.). 
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1 INTRODUCTION 



5-meson decays proceed dominantly through b — > c transitions, while the tree-level b — > it and 
one-loop 6 — ► s transitions are suppressed. In the Standard Model (SM), the branching fractions 
for b — > u and 6 — > s transitions are approximately 1% — 2% [1, 2, 3]. It has been suggested that 
loop transitions are a window on the effects of new physics, as virtual non-SM particles in the loop 
can couple to the quarks [4, 5]. The branching fraction for b — > sg* (g* = gluon) decays could be 
as large as 10% in certain models [4, 5]. 

In recent years, exclusive B decays dominated by b — > sg* (b — > sgg, q = u,d, s) have been used 
to measure the CKM unitarity triangle angle /3. The amplitude S of the sine component of the 
time-dependent CP asymmetry in these decay modes is measured to be systematically shifted low 
relative to the expected values from SM calculations (S ~ sin 2/?), although this shift is currently 
not statistically significant [6]. 

A good understanding of the dynamics of b — > s transitions is needed to make accurate pre- 
dictions of related quantities within the framework of the SM [7]. This understanding currently 
comes from branching fractions and CP measurements of b — > s dominated exclusive decays and 
from inclusive and exclusive b — > 57 transitions. The measurement of inclusive b — ► sg* decays 
would provide additional information to the current picture [8] , and could help us understand the 
discrepancies seen in the measurements of sin 2(3 [7]. 

Previous experimental attempts to measure inclusive b — > sg* decays have been statistically 
limited [9, 10, 11, 12]. The 5-factory experiments present new opportunities to make a significant 
measurement of this process. 

In this paper, we present a preliminary result from a search for inclusive charmless B — > K + X 
and B — > K°X decays, which can in principle be related to the b — > sg* rate. The neutral kaon in 
B — > is reconstructed through the decay -Kg — > 7r + 7r~. We define as signal S — ► i^TX all the 

charmless decays that contain at least one kaon. These decays can occur via b — > s (dominant), 
b — ► it, and 6 — > d transitions. The signal yields are extracted with an unbinned maximum likelihood 
(ML) fit to samples of B — > decays recoiling against fully reconstructed hadronic S decays. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector [13] at the PEP-II asymmetric- 
energy e + e~ collider located at the Stanford Linear Accelerator Center (SLAC). The analysis uses 
an integrated luminosity of 288.5 fb" 1 recorded at the T(4S) resonance (center-of-mass energy 
^s = 10.58 GeV). 

In the BABAR detector, charged particles are detected and their momenta measured by a com- 
bination of a vertex tracker consisting of five layers of double-sided silicon microstrip detectors and 
a 40-layer central drift chamber, both operating in the 1.5-T magnetic field of a superconducting 
solenoid. We identify photons and electrons using a CsI(Tl) electromagnetic calorimeter (EMC). 
Charged particle identification (PID) is provided by an internally reflecting ring imaging Cherenkov 
detector (DIRC) covering the central region of the detector, the average energy loss (dE/dx) in 
the tracking devices and by the EMC. Additional information that we use to identify and reject 
electrons and muons is provided by the EMC and the detectors of the solenoid flux return (IFR). 
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3 ANALYSIS METHOD AND EVENT SELECTION 



The experimental identification of the inclusive b — > sg* decay is complicated by the fact that the 
gluon is a virtual intermediate state with no good experimental signature. We instead rely on 
the hadronization of the primary strange quark into a charged or neutral kaon to identify decays 
dominated by b — > s transitions. The analysis is therefore effectively a search for the decays 
B — > K + X and B — > The momentum p*(K) of the primary kaon in the i? rest frame 

is limited for b — ► c background by the Z)-meson mass, and cannot be larger than ~ 2.3 GeV, 
while p*{K) can be as large as ~ 2.6 GeV for signal B — > i^TX decays. We use this difference as 
the primary signature to look for charmless inclusive B — > i^TV decays, as was first suggested in 
Ref. [14]. 

In this analysis we reject the large e + e~ qq (q = u,d,s,c) continuum background very 
efficiently by reconstructing one of the two B mesons (-B r cco) in e + e~ — > T(4S) — > and 
searching for the S — > ifl signal (B s i gnSL i) recoiling against B Icco . We select a large sample of 

events containing a -B rcco meson which decays into a hadronic final state as B reco — > D V , 
and is fully reconstructed. The system consists of a combination of hadrons containing one, 
three, or five charged pions or kaons, up to two neutral pions, and at most two — > 7r + 7r~ 
candidates. We reconstruct D*~ -> L>°7r~; Z)*° -»■ U°7r ; Z>° -> if+7r~, if + 7r~7r°, J-C + 7r~7r~7r + , 
K°vr+7r~; and D~ -»• K+Tr-TT - , if +7T _ 7r-7r , i^TT - , i^7r _ 7r°, K°7r'vr-7r + . The 5 rcC o candidates 
are characterized kinematically by the energy-substituted mass ities = {\s — p 2 B )z and energy 
difference AE = Eb — \yfs, where (Eb,Pb) is the i?-meson 4-momentum vector, and all values 
are expressed in the T(4S) frame. We require the value of AE to be consistent with zero within 
three standard deviations (a), as measured for each decay mode (10 to 35MeV). We require 
5.25 < mEs < 5.29 GeV, and use this variable in the ML fit described below. We define the purity 
for each -B reC o decay mode as the ratio S/(S + B) measured in control samples, where S is the 
number of reconstructed signal -B reC o and B is the number of background events. We require the 
purity of the selected -B reC o candidates to be at least 20%. In events containing more than one i? reC o 
candidate, we select the decay mode with the highest purity. 

To further reject qq continuum background, we make use of the angle #t between the thrust 
axis of the -B reC o candidate in the T(4S) frame and that of the rest of the charged tracks and neutral 
clusters in the event. The distribution of | cos6*t| is sharply peaked near 1 for combinations drawn 
from jet-like qq pairs, and nearly uniform for the almost isotropic S-meson decays; we require 
|cos#t| < 0.9. Further discrimination from continuum in the ML fit is obtained from a Fisher 
discriminant J 7 , which is an optimized linear combination of four variables: the angles with respect 
to the beam axis of the B momentum and B thrust axis (in the T(4S) frame), and the zeroth and 
second angular moments Lo,2 of the energy flow about the B thrust axis. The moments are defined 
by Lj = J2iPi x l cos Qi\ 3 ! where B% is the angle with respect to the B thrust axis of track or neutral 
cluster i, pi is its momentum, and the sum excludes the B reco candidate. 

We select 2.99 x 10 6 i? reC o candidates with the above criteria, and apply an unbinned ML 
fit to the tit-es and T variables to separate BB events from qq continuum background. We find 
N(B rcco ) = (1.78 ±0.09) x 10 6 , where the uncertainty includes a conservative preliminary estimate 
of the systematics due to the modeling of the data. 

We search for the signal B meson (-B s i gna i) using the charged tracks and the neutral clusters 
that are not part of the B ieco candidate. We reject -B s i g nai candidates containing charged tracks 
compatible with an electron or muon hypothesis, or that contain a reconstructed D meson candidate 
with mass within 30 MeV of the nominal mass. We also require mEs(-Bsignai) > 5.1 GeV. 
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The measured T(45) and -B reC o 4-momentum vectors are used to determine accurately the 4- 
momentum vector of -B s ignab independently of its decay products. We then select charged K + and 
neutral K° s — > ir + ir~ candidates with momentum p*(K) larger than 2.34 GeV, calculated in the 
-Bsignai rest frame. 

The DIRC Cherenkov angle 9 C for charged kaon candidates must satisfy — 5a c < 9 C < +2<r c , 
where a c is the resolution on 6> c , and the upper limit is designed to reject contamination from 
charged pions. To exclude secondary kaons, the distance of closest approach of the K + candidate 
must be within three standard deviations of the -B s ignai vertex, as determined inclusively from all 
tracks recoiling against B rcco . 

We require the mass of Kg candidates to be within ±4.5a of the nominal mass (0.486 < m K o < 

0.510 GeV). The reconstructed K® — > ir + ir~ must have a vertex x 2 probability 7-Vertex > 0.001, and 
its lifetime significance (r/cr r ) must be larger than 3. 

The above selection has an efficiency of (16.1 ± 1.5)% for the decays B — ► K + X with p*(K + ) > 
2.34 GeV, and (6.7 ± 1.1)% for the decays B -► K°X with p*(K°) > 2.34 GeV and reconstructed 
as B -> K° S X, K Q S -► TT+vr". The selected samples contain 246 B -► K + X and 76 B -> K° S X 
candidates. We estimate from Monte Carlo simulation that 10—20% of the selected candidates come 
from b — > c decays, which can produce kaons of momentum higher than 2.34 GeV in the decay of D 
mesons. Contamination from unflavored b —> u and b — > d transitions is negligible for B — > K$X, 
and estimated to contribute to 2.4% of the B — > K + X sample via K/tt mis-identification. 



4 MAXIMUM LIKELIHOOD FIT 



We obtain yields for each decay from an extended unbinned ML fit with the following input ob- 
servables: iues, and p*(K). As described below, the fit is first applied with several probability 
density function (PDF) parameters floating to samples obtained with p*(K) > 1.8 GeV. The signal 
yield is then extracted from a fit to the p*(K) > 2.34 GeV samples, in which the BB background 
yield and p*{K) PDF are fixed to the results of the first fit. 

For each event i and hypothesis j (signal B — ► KX, BB background, continuum background), 
we define the probability density function (PDF) 



K=P J -(m E s i )^(^ i )^(p* 



The likelihood function is 



N 



£ = exp (-]>>,■) J] 



(i) 



(2) 



where Yj is the yield of events of hypothesis j, to be found by maximizing C. N is the number of 
events in the sample. 

The ttt-es and T variables discriminate between BB and qq continuum events. For these vari- 
ables, the same PDF is used for the B — > KX signal and BB background components. The 771es 
PDF for qq continuum is parametrized by an empirical phase-space function [15] of the form 



f(x) oc x \J 1 — x 2 exp — £(1 — x 2 ) 



(3) 



where x = 2mEs/\/*> an d £ is a parameter determined by the fit. For B decays, ttjes is modeled 
by the sum of two Gaussians and the function of Eq. 3 with a different value of £. The T PDF is 
parametrized as a bifurcated Gaussian plus a Gaussian for BB events, and as two Gaussians for 
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qq continuum. The p*(K) PDF is denned over the wider range p*(K) > 1.8 GeV. For the signal 
B — > KX, the PDF is the sum of a phase-space function given by Eq. 3, with x = p*(K) /2.62 GeV, 
and a Gaussian to account for the contribution from exclusive 2-body decays such as B — ► r/'K. 
The parameters of the signal p*{K) PDF will be varied in the evaluation of the systematic errors, 
as the p*{K) spectrum is not well known. The BB background PDF is the sum of three Gaussians, 
two of them used to model the B — > DK and B — ► D*K contributions. The qq component is 
described by the sum of an exponential and a Gaussian. All the PDF distributions are illustrated 
in Figs. 1 and 2. 

The PDF for each variable and each component is initially determined from Monte Carlo (MC) 
simulation. A preliminary ML fit with several free PDF parameters is applied to the sample 
obtained with the relaxed requirement p*{K) > 1.8 GeV. This range of p*(K) includes too much 
background for an accurate determination of the signal yield, but it allows the measurement of the 
yield and the p*(K) PDF for the BB component. The free parameters of this fit are the three 
yields, the fractions of neutral B reco candidates for each component, the size of the p*(K) secondary 
Gaussians for the signal and BB components, the width of the p*(K) main Gaussian for the BB 
component, and the mES exponent parameter for the qq component. The results are illustrated 
in Fig. 1, which shows the projections onto p*(K + ) and p*(Kg) of subsamples enriched with a 
threshold requirement on the signal likelihood computed without the variable plotted. 




BaBar 

Preliminary 



It ft I 



2.4 2.6 2.8 

p*(K + ) (GeV) 



2.4 2.8 

P*(k3 (GeV) 



Figure 1: Projection plots for the p*(K + ) (left) and p*(K$) (right) variables from the fits to the 
p*(K) > 1.8 GeV samples. The projections are obtained with a cut on the signal likelihood (see 
text) retaining about 80% of the signal events. The points are from the data, the full line shows the 
full fit, the dotted line the signal, the short-dashed line the BB background, and the long-dashed 
line the qq continuum background. 



The PDFs determined in the first fit are then used in a second ML fit to the sample obtained 
with p*[K) > 2.34 GeV. Free parameters are the signal and qq continuum background yields, while 
the BB yield is fixed to the fraction of the value measured in the first fit (jp*(K) > 1.8 GeV). 
Systematic errors account for the uncertainties in the fixed BB yield, as determined in the first fit 
and which include the affect of correlations with the signal component. 

Monte Carlo simulated experiments are used to validate the fit procedure, and to evaluate 
possible biases in the yields due to our neglect of small residual correlations among discriminating 
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variables. The bias is determined by fitting ensembles of simulated qq experiments drawn from 
the PDF into which we have embedded the expected number of signal and BB background events, 
randomly extracted from the fully simulated MC samples. The measured biases are listed in Table 1. 

5 RESULTS 

The partial branching fractions are calculated as 

B(B - KX,p*(K) > 2.34GeV) = ^ "^ (4) 

where Ykx is the measured yield, I5 is the fit bias, and e is the reconstruction efficiency The 
results of the fits to the p*(K) > 2.34 GeV samples and the quantities used in the determination 
of the branching fractions are presented in Table 1. The significance is taken as the square root 
of the difference between the value of — 21n£ (with additive systematic uncertainties included) for 
zero signal and the value at its minimum. In Fig. 2 we show projections onto ttt-eS) F and p*(K) of 

Table 1: Number of events to fit A?" can( j, fitted signal yield Ykx m events (ev.), measured bias 
Yj, (see text), detection efficiency e, significance S (with systematic uncertainties included), and 
measured partial branching fraction B for each mode. The first errors are statistical and the second 
are systematic. The quantity in parentheses is the 90% C.L. upper limit for the branching fraction 
B(B -► K°X). 



Mode N cand Y KX (ev.) Y b (ev.) e (%) S (a) B (10" 6 ) 



B 


K+X 


246 


58.4+i° 7 5 


2.2 


16.1 


6.0 


IOR+37+31 
J-yO-34-30 


B 


-^K°X 


76 


21.11$ 


2.8 


6.7 


3.1 


154^+55 (< 266 ) 



subsamples enriched with a threshold requirement on the signal likelihood computed without the 
variable plotted. 

6 SYSTEMATIC STUDIES 

We determine systematic uncertainties affecting the measurement of the yields, the estimation of 
the selection efficiencies, and the measurement of the number of -B rcC o candidates. The systematic 
uncertainties are summarized in Table 2. 

The signal yield systematic errors arise from the fixed BB yields (31.0 ± 6.8 events for the 
charged decay mode and 13.8 ± 4.2 events for the neutral mode), which are varied within their 
uncertainties; the fit bias correction, for which we assign as systematic error the quadratic sum 
of the statistical uncertainty on the correction and one half of the correction itself; and the PDF 
parameter uncertainties, which are left free one by one in the fit and the variation in signal yield 
recorded. The dominant contribution to the PDF parameter uncertainties arises from the poorly 
known signal p*(K) spectrum. We evaluate this uncertainty by floating in the fit the relative size 
and the mean of the Gaussian used in the description of the signal p*(K) PDF. 
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Figure 2: Projection plots for tties (left), T (center), and p*(K) (right) from the fits to the p*{K) > 
2.34 GeV samples. The top plots are for the B — ► K + X decay, and the bottom plots for B — > KgX. 
The projections are obtained with a cut on the signal likelihood (see text) retaining about 85% 
of the B — > ET + X signal events and 75% of the 5 — ► -Kg-X" signal events. The points are from 
the data, the full line shows the full fit, the dotted line the signal, the short-dashed line the BB 
background, and the long-dashed line the qq continuum background. 

Uncertainties on the selection efficiencies are dominated by the statistics of the inclusive B — > 
KX Monte Carlo samples. We also include 0.5% uncertainty per track, 2.1% for the Kg, and 2.4% 
for the K + particle identification criteria. 

The uncertainty in the number of fitted B Teco candidates is taken from the results of that 
fit (5%). 

7 CONCLUSIONS 

We have presented preliminary results for a study of inclusive charmless B — ► K + X and B — ► K°X 
decays, recoiling against fully reconstructed hadronic B decays from T(4S) decays. We measure 
the partial branching fractions for charged and neutral kaons with momentum above the end-point 
for b —> c backgrounds (p*(K) > 2.34 GeV): 

B{B -> K+X,p* > 2.34 GeV) = (196i^(stat.)±^(syst.)) x 1(T 6 , and 

B(B -> K°X,p* > 2.34 GeV) = (154t^(stat.)i^(syst.)) x lO" 6 (< 266 x 1(T 6 at 90% C.L.). 
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Table 2: Systematic uncertainties for the B — > If + A and i? — ► If °A decay modes. The multiplica- 
tive errors are fractional, and apply to the efficiency and to the -B reC o counting, while the additive 
errors are in units of events and apply to the signal yields. 





B -► K+X 


B -► 


Multiplicative errors \/o) 






MP pff 


Q 4 


1 fi 1 

1U. ± 


Tracking eff/qual 


0.5 


1.0 


K° s eff 




2.1 


Kaon FID 


z.4 




Number B reco 


5.0 


5.0 


Total multiplicative (%) 


10.9 


17.0 


Additive errors (events) 






Fixed b — > c yield 


+6.0 
-5.6 


+4.1 
-3.5 


PDF parametrization 


+2.6 
-2.4 


+3.9 
-0.8 


Fit bias 


±1.2 


±1.4 


Total additive (events) 


+6.6 
-6.2 


+5.8 
-3.9 



Known exclusive charmless two-body decays, dominated by the decays B + — ► rj'K + and 5° — > 
rj'K , account for approximately 60% of these branching fractions. Similar two-body decays with 
a K* meson and three-body decays, such as B + — > K + K~K + and B° — > K + K~K°, probably 
account for much of the remainder. 

A theoretical model is necessary to extrapolate these results to the full p*(K) spectrum, and 
ultimately to extract a measurement of the inclusive b — ► sg* branching fraction. Completing this 
extrapolation and assigning the theoretical systematic error from the shape of the p*(K) spectrum 
is the focus of ongoing effort. 
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